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Abstract 

We microscopically analyze thermal spin pumping mediated by magnons, 
at the interface between a ferromagnetic insulator and a non-magnetic 
metal, in the semiclassical regime. The generation of a spin current is dis- 
cussed by calculating the thermal spin transfer torque, which breaks the 
spin conservation law for conduction electrons and operates the coherent 
magnon state. Inhomogeneous thermal fluctuations between conduction 
electrons and magnons induce a net spin current, which is pumped into 
the adjacent non-magnetic metal. The pumped spin current is propor- 
tional to the temperature difference. When the effective temperature of 
magnons is lower than that of conduction electrons, localized spins lose 
spin angular momentum by emitting magnons and conduction electrons 
flip from down to up by absorbing all the emitted momentum, and vice 
versa. Magnons at the zero mode cannot contribute to thermal spin pump- 
ing because they are eliminated by the spin-flip condition. Consequently 
thermal spin pumping does not cost any kinds of applied magnetic fields. 
We have discussed the distinction from the theory proposed by Xiao et 
al. [Phys. Rev. B, 81 (2010) 214418], Adachi et al. JPhys. Rev. B, 83 
(2011) 094410], and Bender et al. [arXiv:1111.2382] . Q 

1 Introduction 

Recently spintronics has developed a new branch of physics called spin caloritron- 
ics [l][2], which combines thermoelectrics with spintronics. Spin caloritronics 
has been attracting a special interest because of potential applications to green 
information and communication technologies [5]. The central theme is the uti- 
lization of thermal fluctuations as well as spin degrees of freedom in order to 
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induce a (pure) spin current. Thus establishing methods for the generation of a 
spin current by using thermal difference, without any kinds of applied magnetic 
fields, is a significant issue. 

In the previous work we have studied quantum spin pumping mediated 
by magnons under a time-dependent transverse magnetic field at the interface 
between a ferromagnetic insulator and a non-magnetic metal. There the fer- 
romagnct act as a source of spin angular momentum; magnon battery named 
after the spin battery [S]. The applied time-dependent transverse magnetic 
field acts as a quantum fluctuation to induce a pumped net spin current under 
a thermal equilibrium condition. Spin angular momentum is exchanged be- 
tween conduction electrons and localized spins via magnons accompanying the 
exchange interaction at the interface. The interface is defined as an effective 
area where the Fermi gas (conduction electrons) and the Bose gas (magnons) 
coexist to interact; the width of the interface is supposed to be of the order of the 
lattice constant [S]. In addition, the pumped net spin current has a resonance 
structure as a function of the angular frequency of the applied transverse field, 
which is useful to enhance the spin pumping effect induced by quantum fluc- 
tuations. Here it should be stressed that magnons accompanying the exchange 
interaction cannot contribute to spin pumping without quantum fluctuations. 
That is, quantum fluctuations (i.e. time-dependent transverse magnetic fields) 
are essential to quantum spin pumping mediated by magnons. 

In this paper, we microscopically propose an alternative mechanism for the 
generation of the spin current without any kinds of applied magnetic fields 
(i.e. quantum fluctuations) ; thermal spin pumping [7] . Inhomogeneous thermal 
fluctuations, i.e. the temperature difference, between conduction electrons and 
magnons induce a net spin current, which is pumped into the adjacent non- 
magnetic metal. This method can be viewed as an alternative way for the local 
spin injection. 

We assume the local equilibrium condition [8]; since the relaxation times in 
the localized spins (i.e. magnons) and conduction electrons subsystems are much 
shorter than the lattice relaxation time [5J [TOl [H] , the reservoirs become ther- 
malized internally before they equilibrate with each other. Therefore we may 
assume that during the relaxation process, conduction electrons and magnons 
can be described by their effective local temperatures; Tg and [71 [12]. Ac- 
cording to Xiao et al. [13], the condition (i.e. temperature difference) can be 
generated by a temperature bias applied over the ferromagnetic film. 

The theoretical setup [12] is almost the same with our previous work [4] 
except the point that applied magnetic fields are not essential; in particular, 
transverse magnetic fields are absent. We consider a ferromagnetic insulator 
and non-magnetic metal junction shown in Fig. [T] where conduction electrons 
couple with localized spins S(x, i), x = {x,y,z) G R'^, at the interface; 



-' xG (interface) 

The exchange coupling constant reads 2J, and the lattice constant of the fer- 
romagnet is oq. In this paper, we take h = 1 for convenience. The magnitude 
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Figure 1: A schematic picture of thermal spin pumping mediated by magnons. 
Spheres represent magnons and those with arrows are conduction electrons. 
When the effective temperature of magnons (Tm) is lower than that of conduc- 
tion electrons (Ts), localized spins lose spin angular momentum by emitting a 
magnon and conduction electrons flip from down to up by absorbing the mo- 
mentum, and vice versa. The interface is defined as an effective area where the 
Fermi gas (conduction electrons) and the Bose gas (magnons) coexist to inter- 
act; J ^ 0. In addition, conduction electrons cannot enter the ferromagnet, 
which is an insulator. 



of the interaction is supposed to be constant and we adopt the continuous limit 
also in the present study. Conduction electron spin variables are represented as 



where are the 2x2 Pauli matrices; [cr-', ct''] = 2iejki<J^ , {j,k,l = x,y,z). 
Operators c^/c are creation/annihilation operators for conduction electrons, 
which satisfy the (fermionic) anticommutation relation; {c^(x, t), c^(x', i)} = 
Snx5{x - x'). 

We focus on the dynamics at the interface where spin angular momentum 
is exchanged between conduction electrons and the ferromagnet. We suppose 
the uniform magnetization and thus localized spin degrees of freedom can be 
mapped into magnon ones via the Holstein-Primakoff transformation; 




(2) 



(cV^c)/2 



(3) 



S+{x,t) = S'^(x,t) +iS'«(x,t) 



(4) 
(5) 
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S- (x, t) = 5'^ (x, t)+iSy (x, t) (6) 
= V2§a^{K,t) + 0{S-^^^), (7) 

5'^(x,i) = ^-af(x,t)a(x,<), (8) 

S* = S/uq, where operators a' /a arc magnon creation/annihilation operators 
satisfying the (bosonic) commutation relation; [a{x,t),a'{x' ,t)] = 5(x — x'). 
Up to the 0(5*) terms, locahzed spins reduce to a free boson system. Conse- 
quently in the quadratic dispersion (i.e. long wavelength) approximation, the 
localized spin with the applied magnetic field along the quantization axis (z-axis) 
is described by the Hamiltonian Hmag, 

^mag= / dxa1'(x,t)(- ^+B)a(x,i), (9) 

Jxe (interface) ^ ^ 

and the Hamiltonian, 7icx{= H^k + T~^'c^)i can be rewritten as 

H^^ = -JSl dx c''{x,t)a'-c{x,t), (10) 

'J x£{intcrfacc) 



n'^^ = -JalJ^ f dx[a''{x,t)c''{x,t)a+c{x,t)+a{x,t)c\x,t)a-c{x,t)]. (11) 

V ^ Jxe(intorfaco) 

The variable m represents the effective mass of a magnon. We have denoted a 
constant applied magnetic field along the quantization axis as B, which includes 
g-factor and Bohr magneton. Let us mention that though we formulate the 
thermal spin pumping theory with B for generalization, in this paper we finally 
take i? = in sec. [3] and discuss the thermal spin pumping effect in sec. SI 
The total Hamiltonian of the system (interface), H, is given as 

n = H,„ag + + Hci, where (12) 



G (interface) 



dx c^(x, t) 



c(x,t). (13) 



The variable md denotes the effective mass of a conduction electron. Eq. (|13p 
shows that 



JS 



(14) 



acts as an effective magnetic field. 

The dynamics at the interface is described by the Hamiltonian Ti'^-^ ; eq. 
(|11[) shows that localized spins at the interface lose spin angular momentum by 
emitting a magnon and a conduction electron flips from down to up by absorbing 
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the spin angular momentum (see Fig. [T]), and vice versa. This Hamiltonian Wg^, 
which describes the interchange of spin angular momentum between localized 
spins and conduction electrons, is essential to spin pumping. Therefore we 
clarify the contribution of magnons accompanying this exchange interaction to 
spin pumping. This is the main purpose of this paper. Here it should be noted 
that we treat localized spins as not classical variables [13] but magnon degrees 
of freedom. As the result, we can microscopically capture the (non-equilibrium) 
spin-flip dynamics on the basis of the rigorous quantum mechanical theory. 

This paper is structured as follows. First, through the Heisenberg equation of 
motion, the thermal spin transfer torque which breaks the spin conservation law 
for conduction electrons is defined in sec. [5] Second, we evaluate it through the 
Schwinger-Keldysh formalism at finite temperature in sec. [3] Last we discuss 
why thermal spin pumping does not cost any applied magnetic fields in sec. IH 
with pointing out the distinction from the farsceing work by Adachi et al. |12| . 



2 Thermal spin transfer torque 
2.1 Definition 

The thermal spin transfer torque (TSTT) [Tl[Tl[Tn], T^"", is defined as the term 
which breaks the spin conservation law for conduction electrons; 

pt + ^-jl^V, (15) 

where the dot denotes the time derivative, js is the spin current density [TB], 
and jOg represents the z-component of the spin density. We here have defined 
the spin density of the system as the expectation value (estimated for the total 
Hamiltonian, "H); 

pl ^ (16) 

In this paper, we focus on the z-componcnt of the TSTT. 

Through the Heisenberg equation of motion, the z-componcnt of the TSTT 
is defined as 

.„, 

This term arises from T-L'^^, which consist of electron spin-fiip operators; 

V^iplM/^- (18) 

Thus, eq. ([T5|) shows that the TSTT (T/ > 0) can be understood as the number 
density of conduction electrons which flip from down to up per a unit of time 
[8], and vice versa. In addition, the TSTT operates the coherent magnon state 

[13. 
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2.2 Pumped net spin current 

In this subsection, we clarify the relation between the TSTT and the pumped 
net spin current. As discussed in the last subsection, the spin conservation law 
for conduction electrons is broken due to the interaction "Hg^; 

Ps+v•j| = 7;^ (19) 

Thus one cannot simply view the time derivative of the spin density for conduc- 
tion electrons, pi, as the spin current density. 

In respect to Planck's constant (we here partially recover fi,), the time deriva- 
tive of the spin density and the TSTT satisfy the relation [T51 [12] ; 

4 = 0{h). (20) 

Therefore is negligible in comparison with at the semiclassical regime, 
where our interest lies. As the result, the spin continuity equation, eq. ([T9|). 
becomes 

V^y-M- (21) 

Then by integrating over the interface, we can evaluate the pumped net spin 
current, / 'f; ■ dSintorface; 

/xe(i„torfacc)rf^V.jf (22) 
interface- 

(23) 

In addition, conduction electrons cannot enter the ferromagnet, which is an 
insulator |20j . Thus the net spin current pumped into the non- magnetic metal 
can be calculated by integrating the TSTT over the interface, eq. ([23]) . 

From now on, we focus on Tg and qualitatively clarify the behavior of the 
thermal spin pumping effect mediated by magnons, at room temperature in the 
semiclassical regime, in sections [3] and S) 



(interface) 



2.2.1 The spin continuity equation for the whole system 

It will be useful to point out that the spin conservation law for localized spins 
(i.e. magnons) is also broken. The magnon continuity equation for localized 
spins [21] reads 

Pfn + V-jS, (24) 

where jm is the magnon current density, and p^-^ represents the z-component of 
the magnon density. We have defined the magnon density of the system also as 
the expectation value (estimated for the total Hamiltonian, T-L); 

pl,^{a^a). (25) 
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In addition, we call 7^ the magnon source term [21] , which breaks the magnon 
conservation law. This term arises also from 'H'^^; 

U^[pl,M/i. (26) 

Through the Heisenberg equation of motion, the magnon source term can 
be determined and it satisfies the relation; 

U = V. (27) 

Then the z-component of the spin continuity equation for the total system (i.e. 
conduction electrons and magnons) becomes 

P?otal + V-jLtal=0, (28) 

where the density of the total spin angular momentum, p^otai' defined as 

Ptotai = pI - Pm, (29) 

and consequently the z-component of the total spin current density, jfojab be- 
comes 

Jtotal = Js ~ Jm (>^0) 

(note that, 5^ = S* — a^a, via the Holstein-Primakoff transformation in sec. [T]). 
The spin continuity equation for the whole system, eq. (|28p , means that though 
each spin conservation law for electrons and magnons is broken (see eqs. (jl9p 
and ([M|). the total spin angular momentum is, of course, conserved [S]. 

2.2.2 The work by Bender et al. 

Last, let us mention a recent preprint [H] by Bender et al., where the authors 
consider a similar problem. We have chosen a different definition of the pumped 
spin current, for reasons now explained. 

Though they have simply recognized the time derivative of the spin density 
for localized spins. 

Pi., (31) 

as the spin current it reads 



pf/'-^pf +V-(jf-jfJ. (32) 

Thus it is clear that even when the total spin angular momentum is conserved 
(eq. (123)), is not directly related to the spin current itself, jf . That is, p^ 
includes other contributions arising from pi and as well as jf . Therefore the 
definition of the pumped spin current by Bender et al. j22j is, in any regime, 
inadequate to their and our case; the mixture of the Bose (magnon) gas and 
Fermi (conduction electron) one. 

That is why, we have adopted different definition of the pumped spin current, 
eq. (ESI, and evaluate the TSTT. 



^ Note that we have adopted our notation. The variable dSf^/dt in Ref. |22| corresponds 
to Pm- 
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3 Schwinger-Keldysh formalism 



The interface is, in general, a weak coupling regime |23j : the exchange inter- 
action, J, is supposed to be smaller than the Fermi energy and the exchange 
interaction among ferromagnets. Thus T-L'^^ can be treated as a perturbative 
term. 

Through the standard procedure of the Schwinger-Keldysh (or non-equilibrium) 
Green's function [H [lH [55] , the Langreth method gZl [Ml US] , the TSTT can 
be evaluated as 

^ [^t^k2,W2^ki,wi^i^ki+k2,wi+W2 ~ ^i%k2,cJ2^ki,wi^i!ki+k2,wi+LJ2] + ^((^) 

The variable Q^^^^ is the fermionic lesser (greater) Green's function, and G^'-^^ 
is the bosonic one. We here have taken the extended time defined on the Keldysh 
contour [26, 27, ,28], c, on the forward path C-..; c = c_>-|-c<_. Even when the time 
is located on the backward path c^ , the result of the calculation does not change 
because each Green's function is not independent; Q'^ — = Q> — where 
gr(a) represents the retarded (advanced) Green's function [21]. This relation 
comes into effect also for the bosonic case 

Each Green's function reads as follows [IS]; 

G<^ = -27nMLo)Siu:~LUu), (34) 
Gl^ = -27r*[l + /B(c^)]5(c^-c^k), (35) 
GaM.ij = 27ri/F(cj)(5(w - w<j,k), (36) 

SaM.u = -2m[l - fF{Lu)]S{LO~LO„,^), (37) 

where the variables (w) and (uj) are the Bose distribution function and the 
Fermi one. The energy dispersion relation reads Wk = Dk^ + B and w^^k = 
FP - {JS + B/2)a- n, where D = l/(2m), F = l/(2mei), a = +l,-l(=t,;), 
and jjL denotes the chemical potential; /i(T) = ep — {nk-QTY / {12e-p) + ©(T^). 
The variable ep represents the Fermi energy. 
Consequently, eg. ([55]) can be rewritten as 

X S{uji - Wki)^(w2 - a;^,k2)5(t^i + - t^i,ki+k2) 

X {[1 + /b(cji)]/f(c^i +CJ2)[1 - fF{L02)] - /b(wi)/f(^2)[1 ^ /f(wi +CJl)]} (38) 

= A7:J^alS J J (^'^Ki + wt,k2 - ^i,ki+k2) 

X [fpii^ki + W^,k2)[l " /F(wf,k2)] + /B(wki)/F(wki + Wt,k2) - /B('^ki)/F('^t,k(^^) 
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Spin-flip condition 

The delta function in eq. p9p represents the condition for spin-flip between 
conduction electrons and magnons. The modes (i.e. ki and which do not 
satisfy this condition cannot contribute to thermal spin pumping. 

The delta function reads 



S{ {D - F)kf - 2Fki ■ k2 - 2JS 



1 



2Ffcifc2 



5 ( cos6' 



{D - F)kl - 2JS 
2Ffcifc2 



(40) 
,(41) 



where cos6' = ki •k2/(fcifc2). Eq. (UD) holds true on the condition; ki ^ 
0, ^2 ^ 0, and F ^ ^. This condition can be justified because the zero-mode 
for conduction electrons (ki = 0) originally cannot contribute to spin pumping 
which is the low energy dynamics; in order to excite the zero-mode so as to 
become relevant to spin pumping, it costs vast energy which amounts to the 
Fermi energy. Such a (relatively high energy) dynamics is out of the system 
we focus on, %. In addition, when the zero mode for magnons (fci = 0) is 
substituted into eq. (j40| . it gives zero because of the finite effective magnetic 
fields JS{y^ 0). Thus the zero- mode of magnons also originally cannot contribute 
to spin pumping and are eliminated. Then we arc allowed to calculate eq. (|39p 
on the condition; ki and ki ^ 0. 

Consequently by using eq. (|4T|) . the TSTT (eq. ([39)) ) can be rewritten as 



iiT^DF^ 



dfci / dkiX^kuki) 
dkiVih), 



NO 

's ' 



(42) 

(43) 
(44) 



where 



1 



(1 - ^)kf ~ 2JS 
2j^hk2 



■ kiki 



^(kl+B)/T^ _ I g(fci-JS-B/2-l+7r2T2/12)/T, ^ i 
1 1 



1 1 



(45) 



We here have defined a variable, C = cos6', and have introduced dimension- 
less variables; ki = y^D/epfci, ki = \/ Fje-pki^ B = B/ep, J = J/ep, 1111(8) = 
Tm(s)/TF = fcBr,ii(s)/eF, where fee denotes the Boltzmann constant. The vari- 
able Tj[i(s) is the effective local temperature of magnons (conduction electrons) 
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[3 HI [13], and 



7:^(^1,^2) (46) 

represents the dimensionless TSTT in the wavenumber space for magnons and 
conduction electrons; 

VCki) (47) 

denotes the dimensionless TSTT in the wavenumber space for magnons, fci, 
after integrating over the wavenumber space for conduction electrons. fc2- Both 
quantities, 7^^(fci,fc2) and 7^^(fci), describe the exchange interaction (J) and 
the temperature (Tjjj(s)) dependence of the TSTT. 

We set each parameter, as a typical follows [131 Uni 130] ; ep = 5.6 

eV, B/eF = 0, F = 4 eV A2, D = 0.3 eV A^, S = 1/2. Here it should be noted 
that we do not apply magnetic fields along the quantization axis; 

B = 0. (48) 




Figure 2: The temperature difference dependence of the dimensionless TSTT, 
Tg^, and the corresponding schematic pictures. Each parameter reads J — 0.002 
and Ts = 300 K. When the effective temperature of magnons is lower than 
that of conduction electrons, localized spins at the interface lose spin angular 
momentum by emitting magnons and conduction electrons flip from down to up 
by absorbing the momentum (a), and vice versa (b). 



4 Thermal spin pumping effect 

Fig. [2] shows that under the thermal equilibrium condition where temperature 
difference does not exist between ferromagnet and non-magnetic metal, spin 
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Figure 3: (a) The spin-flip condition via magnons; z(fci,fc2) = C' = 
[{1-F/D)kl-2JS]{2^/F/Dkik2)~^, where J = 0.002. Magnons at (near) 
the zero-mode cannot contribute to thermal spin pumping because they do not 
satisfy the spin- flip condition, eq. (|40| . (b) The TSTT in the wavenumber 
space for conduction electrons and magnons, T/{ki, k2). Each parameter reads 
J = 0.002, Ts = 300 K, and - = 1.2 K. A sharp peak exists on the Fermi 
wavenumber. (c) The TSTT in the wavenumber space for magnons, l'^{ki); the 
condition is the same with (b). The higher the effective magnon temperature 
becomes, the longer wavenumber of magnons becomes relevant to thermal spin 
pumping. 



currents cannot be pumped because of the balance between thermal fluctuations 
in ferromagnet and those in non- magnetic metal [71 112| I13| . In addition, it can 
be concluded that the pumped spin current is proportional to the temperature 
difference between the magnon and conduction electron temperatures (i.e. Tg — 
Tin); when the effective temperature of magnons is lower than that of conduction 
electrons (see Fig. [2] (a)), localized spins at the interface lose spin angular 
momentum by emitting magnons and conduction electrons flip from down to 
up by absorbing all the emitted momentum jSj, and vice versa (see Fig. [2] (b)). 
This result exhibits the good agreement with the work by Xiao ct al. [T3] ; they 
have reached this result by combining the spin pumping theory proposed by 
Tscrkovnyak ct al. [31] with the Landau-Lifshitz-Gilbcrt equation. 

Figs. [3] (a) and (c) show that magnons at (near) the zero-mode cannot 
contribute to thermal spin pumping because they do not satisfy the spin-flip 
condition between conduction electrons and magnons, due to the finite effective 
magnetic field JS. (see eqs. (gOj), (|45|), and Fig. [3] (a) ). 

The distinction from the work by Xiao et al. and Adachi et al. 

Let us mention that we have set B = Q. That is, a spin current can be 
generated via the thermal spin pumping effect without any applied magnetic 
fields. This point cannot be obtained by Xiao et al. |13| . The pumped spin 
current is proportional to the temperature difference between the magnon and 
conduction electron temperatures; inhomogcneous thermal fluctuations induce a 
net spin current [7] . This is the main difference from the quantum spin pumping 
effect [4]. 

Last we should discuss the distinction from the important work by Adachi et 
al. [12j , with emphasizing that they have already studied thermal spin pumping 
via magnons before our study. They have pointed out that the approach by using 
the stochastic Landau-Lifshitz-Gilbert equation coupled with the Bloch equation 
is equivalent to the one by the Schwinger-Keldysh formalism (i.e. linear-response 
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theory) in the classical regime where quantura fluctuations are negligible. This 
fact has already been confirmed also by the numerical calculation [32] • Though 
they have studied the thermal spin pumping effect mediated by magnons via 
the Schwingcr-Keldysh formalism by the same procedure with our work, un- 
fortunately we have doubts the validity of their calculation; with reflecting the 
statistical properties, the Keldysh Green's function [55] for fermions {G^) should 
be[25] 

^k<. + GuJ = 2* Im^^_„tanh(/3w/2), (49) 

not 2i ImQl^ ^coth(/3aj/2) [T^- The variable Q^^ ^ denotes the fermionic retarded 
Green's function and f3 is defined as /3 = l/{kBT). That is, the fermionic 
Keldysh Green's function is different from the bosonic one. In addition, we 
would like to mention that though they have taken a classical approximation, 
we have discussed the thermal spin pumping effect in the scmiclassical regime. 
Moreover, we stress that thermal spin pumping does not cost any applied mag- 
netic field, magnetic fields along the quantization axis nor transverse magnetic 
fields, because magnons at the zero-mode are eliminated because of the spin-flip 
condition, eq. PO]) . 

5 Summary and discussion 

We have qualitatively studied thermal spin pumping mediated by magnons in 
the scmiclassical regime. Pumped spin currents are proportional to the temper- 
ature difference between conduction electrons and magnons. That is, inhomo- 
geneous thermal fluctuations induce a net spin current; when the effective tem- 
perature of magnons is lower than that of conduction electrons, localized spins 
lose spin angular momentum by emitting magnons and conduction electrons flip 
from down to up by absorbing the momentum, and vice versa. Thermal spin 
pumping has the advantage that it docs not cost any kinds of applied magnetic 
fields because magnons at the zero mode are eliminated due to the spin-flip con- 
dition. This fact will be useful for potential applications to green information 
and communication technologies; spin currents can avoid Joule heating. 

Though the behavior of the thermal spin pumping effect mediated by magnons 
can be qualitatively captured by calculating the TSTT, we recognize that the 
theoretical estimation for the width of the interface, so called proximity effects, 
is essential for the quantitative understanding. In addition, we are also inter- 
ested in the contribution of phonons and that of magnons under a spatially 
nonuniform magnetization to spin pumping. 
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